Version of 4 October 2012 



Energy separation in oscillatory Hamiltonian 
systems without any non-resonance condition 

Ludwig Gauckler 1 , Ernst Hairer 2 , Christian Lubich 3 

1 Institut fur Mathematik, TU Berlin, Strafie des 17. Juni 136, D-10623 Berlin, Germany. 
E-mail: gaucklcr@math.tu-berlin.de 

2 Section de mathematiques, 2-4 rue du Lievre, Universite de Geneve, CH-1211 Geneve 4, 
Switzerland. E-mail: Ernst. Hairer@unige.ch 

3 Mathcmatischcs Institut, Universitat Tubingen, Auf dcr Morgcnstcllc, D-72076 Tubingen, 
Germany. E-mail: Lubich@na.uni-tuebingen.de 



Abstract: We consider multiscale Hamiltonian systems in which harmonic os- 
cillators with several high frequencies are coupled to a slow system. It is shown 
that the oscillatory energy is nearly preserved over long times e~ N for arbitrary 
N > 1, where e _1 is the size of the smallest high frequency. The result is uniform 
in the frequencies and does not require non-resonance conditions. 



1. Introduction and statement of the main result 

As has been discussed, e.g., in [2,3,4,5,15], the problem of slow rates of ther- 
malization in statistical mechanics can be understood by studying multiscale 
problems of the following type: For momenta p = (po,pi, ■ ■ ■ ,p n ) and positions 
Q = (?0) 9ij • • • > Qn) with pj, qj e K dj consider a Hamiltonian that couples high- 
frequency harmonic oscillators with a Hamiltonian of slow motion, 

H(p,q) = H M (p,q) + # s i ow (p, q), (1) 

where the oscillatory and slow-motion energies are given by 

n 

^(P-ql^^lftP+^lftl 2 ), if slow (p,q) = ibo| 2 + E/(q) (2) 

i=i 

with high angular frequencies 

LJj>-, 0<e<l. (3) 

The coupling potential U (q) is assumed smooth with derivatives bounded inde- 
pendently of the small parameter e. 
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Weakly anharmonic systems similar to (l)-(3) have been studied in numerous 
publications, from different viewpoints and with different mathematical tech- 
niques, ever since the seminal paper by Fermi, Pasta & Ulam ([9], reprinted in 
[fO]). For some recent work, see [1,8,10,12,14] and references therein. The sys- 
tem (l)-(3) has been considered in [4] with the double motivation of explaining 
slow rates of thcrmalization and the realization of holonomic constraints. 

The equations of motion are 



where Vj denotes the partial derivative with respect to qj and luq = 0. We are 
interested in the energy exchange (or lack thereof) between the fast and the 
slow system, where in contrast to previous work we assume no condition on the 
frequencies other than the lower bound (3). 

Theorem 1. Fix an energy bound E > 0, an integer N > 1, and further a radius 
p > and a set K C M. d ° such that the potential U has bounded derivatives of 
all orders in a p -neighbourhood of K x x ■ • • x 0. Then there exist C > and 
e* > such that the following holds for < e < e* : Whenever the frequencies 
satisfy loj > for j = 1, . . . , n and the initial values (p(0), q(0)) are such that 



then, along the solution of (4) to these initial values, the oscillatory energy de- 
viates from its starting value by no more than 



H„(p(t)Mt)) -^(p(0),q(0))| < Ce 3 / 4 for < i < e~ N , (6) 



provided that qo(t) stays in the set K for such long times. The threshold e* and 
the constant C depend on n and N , on the energy bound E and on bounds of 
derivatives of the potential U . 

We emphasize that no non-resonance condition is required in this result. The 
estimate is uniform in the frequencies u)j > e . This sets Theorem 1 apart from 
existing related results in Hamiltonian perturbation theory. 

For fixed strongly non-resonant frequencies as well as for exactly resonant 
frequencies that satisfy a diophantine non-resonance condition outside a res- 
onance module, the long-time estimate (6) is a particular case of the energy 
separation results shown by Benettin, Galgani and Giorgilli [4]. The proof of 
Theorem 1 is based on modulated Fourier expansions with suitably modified 
frequencies and does not use the canonical transformation techniques of Hamil- 
tonian perturbation theory on which the proofs in [4] are based. The idea of 
modifying frequencies to avoid almost-resonant situations has also been used in 
[2] in a different context where all frequencies are sufficiently close to a single 
high frequency 

The results in [4] are proved over times that are exponentially long in negative 
fractional powers of e in the case of a real-analytic potential. We expect that then 
also Theorem 1 holds over exponentially long times, but we will not investigate 
this question in this paper. See, however, [6] for a corresponding result over 
exponentially long times for a single high frequency proved via modulated Fourier 
expansions. 



j = 0, . . . , n, 



(4) 




(5) 
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We remark that the order e 3//4 can be improved, with the same proof, to e 1- * 5 
for arbitrary S > 0, with e* and C depending additionally on S and deteriorating 
as S — > 0. 

The energy bound if w (p, q) < Const, implies qj = 0(e) for j = l,..,,n. 
Theorem 1 and conservation of the total energy therefore yield that also the 
energy ^\po\ 2 + U(qo, 0, . . . , 0) of the isolated slow system is nearly preserved 
over time e~ N along solutions of (4). 

We present a numerical example in the next section and prove Theorem 1 in 
the remainder of the paper. 

2. Numerical illustration 

We consider the problem with n = 7 frequencies 

euj = 1, 1 + e 2 , 1 + e, 1 + e 3/i , 1 + e 2/3 , 1 + e 1/2 , 2 
for j = 1, ... ,n, and with the potential 

U(q) = - q 2 + (a q + qi + q 2 + 2 q 3 + 3 q 4 + q 5 + <7e + 3 q 7 ) 3 
depending on a parameter a. The initial values are chosen as 

q(0) = (l,0.3£,0.4e,0.7e,-l.l£,0.4e,-0.6e,-0.7e) 
p(0) = (-0.2,0.6,0.7,-0.9,-0.9,0.4,-1.1,0.8). 

The numerical solution is computed with a trigonometric integrator (symplectic 
method of Dcuflhard, see [13, Chap. XIII]), which is applied with constant step 
size At = 0.01 e. Taking the step size twice as large gives nearly identical figures. 




30000 60000 90000 

Fig. 1. Individual oscillatory energies and shifted total oscillatory energy with parameter 
a = 0.5 and e = 0.005. 



We consider two different values for the parameter a. Figure 1 shows the 
individual oscillatory energies Ej = ^(\pj\ 2 +w 2 |g :) | 2 ) as well as the total oscil- 
latory energy on an interval of length 100 000 for a — 0.5 (for a = 1 the 
energies explode) and for e = 0.005. The initial oscillatory energies are close to 
0.22,0.32,0.65,1.05,0.17,0.82,1.3. Wc observe that the oscillatory energies for 
the modes j = 4, 5, 6 remain nearly constant, whereas for the other modes, which 
are closer to a 1:1 and 1:2 resonance, there is a significant energy exchange. In 
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the figure the constant value 2.3 is subtracted from the total oscillatory energy, 
so that it fits better into the picture. We repeated the experiment with other 
values of e and found that the maximal deviation (6) of the total oscillatory 
energy is 4.56 ■ 10" 1 for s = 0.02, 1 .85 • 10" 1 for e = 0.01, and 9.54 • f0~ 2 for 
e = 0.005, which indicates an 0(e) -behaviour. 
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Fig. 2. Individual oscillatory energies and total oscillatory energy with parameter a = e. 
The upper picture is for e = 0.02, the middle picture for e = 0.01, and the lower picture for 
e = 0.005. 



Figure 2 shows the same experiment for the parameter a = e and for three 
different values of e. Here, the maximal deviation (6) of the total oscillatory 
energy is 3.95 • 10~ 2 for e = 0.02, 1.41 • 10~ 2 for e = 0.01, and 4.81 • 10~ 3 for 
e = 0.005. 



3. Modulated Fourier expansion with modified frequencies 



3.1. Changing almost-resonant frequencies to exactly resonant frequencies. We 
deal with almost-resonant frequencies by introducing modified frequencies close 
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to the original frequencies such that almost-resonances among the original fre- 
quencies become exact resonances among the modified frequencies. For a multi- 
index k = k n ) G Z n , we denote ||k|| = fci| + . . . + \k n \. The vector of 
frequencies is a> = (usi, . . . , w„), and we write k ■ uj — Y^j—i kjUij. 

In order to decide whether a linear combination k ■ u; of frequencies is consid- 
ered as almost-resonant or as non-resonant we use the fact that there is a gap 
in the linear combinations of frequencies. 

Lemma 1. There exist /i > depending only on n and N and a G [/i, j] depend- 
ing on the frequencies such that there is no multi-index li G Z™ with ||k|| < N + l 
and£- a < |k-w| < s- a -v. 

Proof. Let us denote the number of all multi-indices k G Z n with ||k|| < N + 1 
by M , which obviously depends only on n and TV. Consider the logarithms to 
base e _1 of the |k • u>| with ||k|| < N + 1, which gives a set {oi, clm}- With 
fj, = 1/(4M + 4), there is an interval of width fj, in [/x, j] that does not contain 
any of the a^. If we denote this interval by [a, a + /x], then there is no k G Z" 
with ||k|| < + 1 and e" Q < |k • w| < e~ a -^. □ 

Let a and fi be as in the previous lemma. We consider multi-indices k G Z n 
with ||k|| < N + 1 and |k • u;| > e _Q_A1 as non-resonant, whereas those with 
|k • u>| < are considered as almost-resonant and collected in the set 

K \k : .." : |k • oj\ < £- a , ||k|| • A • l| . 

This set generates the module M. that contains all finite linear combinations of 
multi-indices in 1Z with integer coefficients. We select a maximum number d < n 
of R- linearly independent vectors k 1 , . . . , k d G 1Z. Their real- linear span contains 
the module M. 

We modify the almost-resonant frequencies ujj to exactly resonant frequencies 
w i = u j 3 = 1) ■■■,n, 

by determining i9 = . . . , as a solution of minimal norm of the (possibly 
underdetermined) system of linear equations 

k i -a; + k i -i? = 0, i = l,...,d. 

We then have 

M<ie- a , (7) 
where 7 depends only on n and N, and 

k-zu = forallkGA^, (8) 

i.e., almost-resonances in the original frequencies become exact resonances in 
the modified frequencies. We further have the lower bound 

|k-w| > if -" for k^M, ||k|| <JV + 1, (9) 

provided that < s < e* , where e* depends only on n and N, i.e., non- resonant 
multi-indices for the original frequencies arc still non-resonant for the modified 
frequencies. It will be convenient to let wq = and ^0 = 0. 
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3.2. Short-time approximation by a modulated Fourier expansion. We let TV be 
a set of representatives of the equivalence classes in Z, n /A4, which are chosen 
such that for each k 6 TV, the norm of k is minimal in the equivalence class [k] = 
k + M, and with k G TV, also -k G TV. We denote JC = {k G TV : ||kj| < TV}. 
We make the approximation ansatz (modulated Fourier expansion) 

9o(t) «X) z i( e ~ a *) ei0t " 7)t ' i = 0,l,...,n. (10) 
keic 

The modulation functions Zj(t) for t = e~ a t with a from Lemma 1 are deter- 
mined such that all their derivatives with respect to r arc bounded independently 
of e and that they yield a small defect 6}f in the modulation equations, which 
are derived by inserting the expansion (10) into the differential equation (4) and 
comparing the coefficients of e 1 ^"^ 7 )*: 

(g7? - (k • zu) 2 ) zf + 2 i(k ■ zu) s~ a ^r + £ ~ 2a -£r ( n ) 

= (2^ - - VJ k U(z) + » k . 

Here V~ k denotes the gradient with respect to z~ k and the modulation potential 
U(z) for z = (zj ) is given as 

JV n 

= c/(z°) + J2 E E 

m=l ji,...,j m =0 k 1 + ...+k'"GA4 

where z° = (z° , 0, . . . , 0) and the indices (ji, k') = (0, 0) are excluded from the 
sum. Note that z*~ are vectors in M. dj and correspondingly 9^ . . . dj m U(z®) is a 
multilinear form on R dj i x • • • x K im ; in contrast VjU denotes the gradient, a 
column vector in M. dj . In the following we let (J) — (0, . . . , 1, . . . , 0) denote the 
jth unit vector in Z n . 

Proposition 1. In the situation of Theorem 1, the solution Qj(t) (j = 0, . . . ,n) 

of (4) admits an expansion 

q 3 {t) = zfie-^e^^ + rM 0<t<e a , (12) 

where the coefficient functions z^(t) satisfy z~ k = z^ and are bounded for < 
T < 1, together with all their derivatives with respect to r up to any fixed order, 
by 

\z^{r)\<C u \zf^(r)\ <duj l <C lS , 
and for all other (j, k), 

^(rJl^dlt^-Ck-w) 3 !- 1 ^". 

The functions z^f satisfy the modulation equations (11) with a defect bounded by 

\Sj(r)\<C 2 s N+1 , 0<r<l. (13) 
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The remainder term is bounded by 

(u>*\ rj (t)\ 2 + Mt)\ 2 f 2 < C,e N+ \ 0<t<e a . 

The constants Ci, 62,63 are independent of e and the frequencies OJj > £ _1 , but 
depend on n and N , on the energy bound E and on bounds of derivatives of the 
potential U . 

This will be proved in Section 4. The proof follows the lines of previous proofs 
for modulated Fourier expansions; compare for example [13, Chap. XIII] and, in 
technically more complicated situations than here, [7,11,12]. 



3.3. The almost-invariant. We introduce the functions y^{t) — z^(e Q f)e^ k ' ro ^, 
which we collect in the vector y = (y k ) for j = 0, . . . , n and k G JC. Since kz& = 
for k e M, the modulation potential satisfies U{y) = U(z). By (11), we have 
with rotated defects d%(t) = $jf(e-°<)e i ( k,OT ) t that 

# + w?!£ = -V7 k W(y) + d£ 1 (14) 

where V~ k denotes the partial derivative with respect to yj . 

Without the defects c? k , equation (14) is again a Hamiltonian system. Since 
the sum in U(y) is over multi-indices k 1 , . . . , k Tn with k 1 + • • • + k m G A4, the 
extended potential IA is invariant under the action of the one-parameter group 
S(6)y = (e^*-^): 

U{S{9)y) = U{y) for all 6 G R, 

because k ■ vj = for k G M. By Nocther's theorem this leads to a conserved 
quantity £ of the system (14) considered without defects. As we now show, this 
quantity is almost conserved in our situation including the defects d k : We have 



- E ' w) yj k Vj k U(y) = - U(S(9)y) 
j=o ke/c 



= 0, 

8=0 



and for 



f(y,y) = -iEE( k -^)y7 k ^ k ( 15 ) 

j=o ke/c 

we thus obtain by multiplying (14) with — i(k ■ vj)y~^ and summing over j and 
k that 

d " 

j t £ (y(t), y(t)) = -i E E ( k ■ ro ) yJ^) # (*)• 
i=o ke/c 

Since the right-hand side is 0(e N+1 ) by the estimates of Proposition 1, £ is an 
almost-invariant. Moreover, £ turns out to be close to the oscillatory energy 
along the solution of (4). 



8 



L. Gauckler, E. Hairer and Ch. Lubich 



Proposition 2. In the situation of Proposition 1 we have for < t < e a 

l%W,y(«))-«(y(o)j(o))| <Ce N+1 

\£(y(t),y(t)) - H u (p(t),d(t))\ < C'e 3 /\ 

The constants C and C are independent of £ and the frequencies cuj > e . 

Proof. The first inequality has already been derived above. Concerning the sec- 
ond inequality, we note that with r = e~ a t we have for k ^ 

d k 

= e^-* (e- a ^r(T) + i(k • ™)zf(r)) = i(k • vj)${t) + 0{e^). 
With the bounds of Proposition 1 it follows that (omitting the argument t) 

eM)=2Y t rf\yW\* + Q^-°) 
i=i 

71 

= 2Y,^\y i / ) \' 2 + 0(e 1 - a ). 
On the other hand, with Proposition 1 wc also obtain 

n 

^( P ,q) = lE(i^f +yj {j) \ 2 + ^\vf ) +vj {j) \ 2 ) +o(e) 

3 = 1 

n 

= 2]>>>f| 2 + 0(e 1 -«). 
i=i 

Since a < |, the result follows. □ 

5.^. Transition from one interval to the next. Propositions 1 and 2 are only 
valid on a short time interval < t < e a . We now consider a second short time 
interval e a < t < 2e a and control in the following proposition the transition in 
the almost-invariant from the first to the second interval. 

Proposition 3. Let the conditions of Proposition 1 be fulfilled. Let z^(e^ a t) 
for < t < e a be the coefficient functions as in Proposition 1 for initial data 
(p(0),q(0)), and let yf(t) = z^(e _ °i)e i ( k - w ) t and y(t) = (y*(i)). Let further 
y(t) = (yj(t)) be the corresponding functions of the modulated Fourier expansion 
for < t < e a to the initial data (p(e Q ), q(e Q )), constructed as in Proposition 1. 
Then, 

\£(y(e a ),y(e a ))-£(y(0)M0))\ <Ce N+ \ 
where C is independent of e and the frequencies uij > e^ 1 . 



The proof will be given in Section 4. 
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3.5. From short to long time intervals. Now we put many short time inter- 
vals together to get the long-time result of Theorem 1. For m = 0,1,2,..., 
let y m {t) contain the summands of the modulated Fourier expansion starting 
from (p(me a ), q(me a )). As long as H u (p(me a ), q(me Q )) < 2E, Proposition 2 
yields for < t < e a 

\£(y m (t),ym(t))-£(y m (0),y m (0))\ < Ce N+1 . 

By Proposition 3, 

\£(y m (e a ),y m (£ a )) -£(y m+ i(0),y m+ i(0))| < Ce N+1 . 

Summing up these estimates over m and applying the triangle inequality yields, 
for < < e a , 

|£(ym(0),ym(0))-£(yo(O),y o (O))| <2(m + i)Ce N+1 . 

By Proposition 2, we have at t = me a + 6 

\£(y m (9),y(0))-HMt)^(t))\<C'e 3 /\ 
Combining these bounds at t and at we obtain for t < s~ N 

|fl«(p(t),q(t))--H'«(p(0),q(0))| < 2Cte N+1 ~ Q + 2C'e^ 4 . 
This yields the bound of Theorem 1. 



4. Proof of Propositions 1 and 3 

4-.1. Construction of the modulation functions. We give an iterative construction 
of the functions Zj (j = 0, . . . ,n and k S JC) such that after M = \(N + 

iteration steps, the defect in equations (11) is of size 0(e N+1 ). We denote by 
z m = ([z^] 1 ™) the mth iterate and distinguish between the following cases: 

1. For j = and k = the first two terms in (11) disappear and we iterate with 
a second order differential equation for [zq]" 1+1 : 

j2r_Olm+l 
-2a u n) V7-O7// m\ 

s — 2 = -V W(z ). 

2. For j ^ and k = ±(j) with the jth unit vector (j) = (0, . . . , 0, 1, 0, . . . , 0), 
the first term in (11) disappears and we iterate using a first order differential 
equation for [zf^] m+1 : 

±2i Wj£ - a ^-^ +e- 2a L j 9 J (16) 

dr dr 2 

= (2w^j - tf)[zf a> ] m+1 - Vf <J> U(z m ). 



We note that by (7) we have for sufficiently small e 



2n7 J e- Q 



< 27. 
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3. In all other cases we iterate with an explicit equation for [z k ] m+1 : 

K - (k ■ -) 2 ) [zf] m+1 + 2 i(k • W ) e-^p- + e- 2a ^§f- (17) 

= (2w j j -&)[zf] m -Vj k U(* m )' 



For j ^ wc note that ±(j) — k ^ .M for k G /C and k 7^ ±(j), and therefore 
we have by (9) 



(k- zo)E~ a 



(k • zaf 



< 2£ A 



-2a 



-(k-w) 



< 2e 



1 — Q + /i 



and 



07? - (k • ^) 2 



< 87^. 



For ,7 = and k ^ jK we have similar estimates directly by (9). 

We need initial values [z£] m+1 (0), £[z%] m+1 (0) and [zf {j) ] m+1 (0) for j ^ 0. 
They are uniquely determined by the equations (for j = 0, . . . , n) 



9i(0) = E #] m+1 (0) , & (0)=£(V Q (0) + i(k • w) [z k ] m+1 (0)) 

(18) 

after inserting (17) for k 7^ ±(j) and replacing for j 7^ the derivative of 
j z ±0>jm+i w jth (16). The starting iterates are chosen as [zo]°( T ) equal to the 
solution at t = e a r of y = -V U(y, 0, . . . , 0) with y(0) = q (0) and y(0) = q (0), 
and for all other (j, k) with k 7^ ±(j), we set [z^]°(t) = 0. The diagonal functions 
j z ±(j)jO are than chosen as constant functions such that (18) is satisfied. 



4-2. Bounds of the modulation functions. We estimate the modulation functions 
on intervals of length r = 0(1), or equivalcntly t = 0(e a ). With the above 
inequalities and the energy bound (5) we obtain by induction on to, for m = 
0,...,M, 

[ Z °r = o(i), [ Z f {j) r = o^j 1 ) 

and for all other (j, k), 

[ 2 *r = o(it^-(k. W ) a |- i e ii k ») ) 

and the same bounds hold for all their derivatives up to any prescribed order. 
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4-3. Bounds of the defect. Wc prove the estimate (13) of the defect in Propo- 
sition 1. The defect [5^} m (r) in inserting the mth iterate z m = ([z^] m ) in the 
modulation equations (11) satisfies 



'j2r Olm+l J2[ Olm 
rrOijri _ ,-2a " l z J " l z 0l 



and for j = 1, 



[^=W]m = - F2in7j . e 



dr dr 



and for all other (j, k) 



rrkim / 2 /i \2\ / r kim+1 r kim 

=-K- -(k-w) ) (J^-] 

In the following we work with the C r ([0, 1]) norm, defined for v = (v^) by 

" d' 
||v||cr = max max ^ £ 



0<r<l 0<Kr /-^ \ dr 1 

" j=0 k£K 



and use the notation 



Az = v = (vf) with < = 2i WjS^zf^ , j = 1, . . . , n 

With this notation and the above formulas for the defect we have 

|[<5j k r| < C||v m+1 - v m || C 2. 

We therefore study v m+1 — v m and show by induction on m that for m < M 
and r = 2M - 2m + 2, 

|| v m+i -v m \\ cr = 0(e mfl ), (19) 

which implies the estimate (13) of the defect in Proposition 1 if we use M > 
(N + l)//z iterations for the construction of the modulation functions. Estimate 
(19) is shown by analysing the iterative construction of modulation functions, 
where the (m + l)th iterates are the dominant terms and the mth iterates the 
non-dominant ones. 

We first consider the diagonal elements [v^^] m . From (18) we obtain after 
inserting (17) for k ^ ±{i) that 

\[vf {j) } m+1 (0) - [vf {j) } m (0)\ < Ce»\\v m - v" 1 - 1 !!^. 

Solving the differential equation (16) by the variation-of-constants formula then 
yields that 



3 



(r) - [vf w ] m (r)| < Ce"\\v m - v'"- 1 ^, < r < 1. 
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Repeated differentiation in (16) further shows that for I < 2M — 2m + 2 and for 
< t < 1 

< Ce^||v m -v m - 1 ]| C €+2. 
Using 2a > p we also bound, for r = 2M — 2m + 2, 

HKP+ 1 - [v°r\\ C r < cs»\\v n - v™- 1 !^. 

Similarly we have from the iteration formulas and the above inequalities that 
|| v m+l _ v m|| Cr < C£/i || v m_ v m-l|| Cr+2 _ 

This proves (19) and the estimate of the defect of Proposition 1. 



4--4- Solution approximation. We consider the Mth iterates of the modulation 
functions (with M = \(N+l)//i\ ) and omit the superscript M on the modulation 
functions and their defects. The truncated modulated Fourier expansion 

^) = ^ k ( £ -^)c i(k - ro)t 
ke/c 

satisfies the perturbed second order differential equation 

qj + ujjqj = -Vji7(q) +dj, j = 0, . . . , n, 

with 

dj(t) = ^ k (e _a i)e i(k - ro)t - VT k ^(z(£- tt t))e i ( k - ro )* + Pj (t) = 0(e N+1 ), 

ke/c keJV\/c 

where p(t) denotes the gradient of the remainder term in the truncated Taylor 
expansion of the potential U about , i®(£~ a i): omitting the arguments t and 
t = e~ a t, and writing q = z° + q, 

N n 

Pi =V i C/(q)-V J -C/(zg)-^ J2 —i - ■•• 9 i W*8)($ </, )• 

m=l j'i,...,j m =0 

After subtracting the differential equation (4), standard estimates using the vari- 
ation of constants formula and the Gronwall inequality then yield the bound of 
Tj(t) = qj{t) — qj{t) for < t < e a as stated in Proposition 1. 



« lv Mih 

nlT 1 [V 3 1 



1 ' air 1 1 
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4.5. Proof of Proposition 3. Let z m = ([z^] m ) and z m = {[z*] m ) be the mth 
iterates of the construction of the modulation functions z and z for initial data 
(p(0),q(0)) and (p(e Q ), q(e Q )), respectively. With the iteration for z m and the 
defect formula for z M the same arguments as in the estimates of the defect yield 
for w m (r) = Az M (l + r) - Az m (r) that 

\\w M \\ c2 =0(e N+1 ). (20) 
The only difference is that the estimate of [iVj' > ] m+1 (0) now reads 

\[wf {j) ] m+1 (0)\ < C^||w m || c2 + Ce N+ \ 

a bound that is obtained by inserting (12) in the defining relation (18) for 
[z^ ] m+1 (0) and by using the estimate of the remainder term from Propo- 
sition 1. Similarly the estimate of w° has to be modified. The bound of Propo- 
sition 3 follows from (20). 
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